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Abstract 

We revisit the question of whether or not one can perform rehable semiclassical QCD 
computations at zero temperature. We study correlation functions with no perturbative 
contributions, and organize the problem by means of the operator product expansion, es- 
tablishing a precise criterion for the validity of a semiclassical calculation. For Nf > N, a 
systematic computation is possible; for Nf < N, it is not. Nf = is a borderline case. 
In our analysis, we see explicitly the exponential suppression of instanton effects at large 
N. As an application, we describe a test of QCD lattice gauge theory computations in the 
chiral limit. 



1 Introduction: Perturbative and Non-Perturbative contribu- 
tions to the OPE 



Over the last few decades, semiclassical computations have proven valuable in understanding 
a variety of effects in quantum field theory and string theory. In four dimensional gauge 
theories, 't Hooft showed in his original work that systematic computations arc possible in a 
Higgs phase[l]. It was soon realized that in QCD at high temperatures, instantons provide 
the leading contribution to ^-dependent quantities, as well as to certain chirality violating 
amplitudes[2, 3]. These methods found applications to the understanding of baryon number 
violation in the Standard Model, the axion potential at high temperatures and the dynamics of 
super symmetric gauge theories. But the question of the utility of instantons at zero temperature 
in real QCD has remained elusive. Even for short distance quantities, such computations arc 
plagued by infrared divergences. These divergences arise precisely in the region of strong 
coupling, suggesting that instantons arc, at best, just one of many important non-perturbative 
effects. Witten, indeed, argued compellingly that at large N, instanton effects are negligible [4]. 

A number of authors in the past have considered possible contributions to physical pro- 
cesses in QCD by studying momentum space Green's functions[5, 6, 7]. Such computations, 
remarkably, are infrared finite, and arguably might represent distinguishable effects. But there 
are reasons to be skeptical. Most strikingly, in the cases which have been studied, these in- 
stanton contributions fall very rapidly with momentum. It is not clear that, e.g. in large 
N, the leading contributions to the corresponding quantities won't swamp the semiclassical 
contributions. 

In this note, we resolve this question by considering a set of correlation functions which 
receive no perturbative contribution but do receive contributions from instantons. They receive 
other non-perturbative contributions as well, as readily seen in the large N limit. By considering 
the structure of the Operator Product Expansion (OPE), we can show that, for particular 
numbers of flavors and colors, instantons provide the most singular contribution to these Green's 
functions at short distances, while for others they are subdominant. In the former cases, as 
we demonstrate self-consistently, a complete semiclassical computation is possible. There is a 
systematic expansion of the most singular part of the coefficient function. The leading instanton 
contribution to the singularity is infrared finite, and corrections can be calculated systematically 
in powers of a^. Dilute instanton gas corrections are infrared divergent, and can be organized 
as corrections to the matrix elements of higher dimension operators, as well as less singular 
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corrections to coefficient functions. 

There is also a satisfying connection to large N arguments. In the large N limit, the 
coordinate space instanton contributions to these correlators are always infrared divergent, and 
the momentum space expression, while formally finite, is exponentially small compared to the 
leading non-perturbative contributions. These latter can be determined systematically in terms 
of (qq). The computation provides an explicit realization of Witten's assertion that instanton 
effects are exponentially suppressed in large N[4]. A similar analysis demonstrates that in e+e~ 
annihilation and deep inelastic scattering, instanton effects at large momentum transfers are 
highly subleading, and cannot be discussed in a systematic fashion (independent of the number 
of flavors and colors). As we will also explain, our observations may provide a test of lattice 
computations which attempt to penetrate the chiral region. 

The rest of this paper is organized as follows. In the next section, we introduce the 
maximally chirality violating Green's functions which interest us. We obtain the condition 
that the instanton contributions are finite and show that, in such cases, there is always a 
power law singularity at short distances. We explain how to organize the computation in 
the language of the operator product expansion. In the cases that the instanton contribution 
is infrared divergent, we show that the infrared divergence corresponds to a correction to the 
matrix elements of operators already appearing in the perturbative OPE, as well as non-singular 
corrections to the unit operator and operators of higher dimension. In the case where the result 
is infrared finite, we demonstrate that the result represents a calculable, singular contribution 
to the coefficient of the unit operator due to instantons. 

In section 3, we discuss corrections and establish the nature of the semiclassical expansion. 
We show that perturbative corrections about the instanton yield corrections controlled by as{x), 
and are systematically small at short distances. The dilute instanton gas corrections are infrared 
divergent, and can again be understood as corrections to matrix elements and non-singular 
corrections to OPE expansion coefficients. In section 4, we consider the large N approximation, 
showing that one can systematically calculate the Green's functions in terms of (qq), and that 
the leading large N contributions in momentum space are exponentially large compared to the 
instanton contribution. In section 5, we consider the case of Nf = N, which is the borderline 
for calculability. Here it is possible for instanton effects to dominate, but only by powers of 
logarithms. We compute the power for two interesting cases. In section 6, we consider the 
application of these observations to e+e^ annihilation and similar processes, demonstrating 
that the instanton effects are subleading, and cannot be organized as part of any systematic 
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semiclassical expansion. In section 7, we discuss applications to a possible test of lattice gauge 
theory computations. Finally, in our concluding section, we outline two directions for further 
work: a study of these questions in CP^ models, and a detailed analysis of the lattice gauge 
theory problem. 

2 Instantons and Maximally Chirality Violating Green's Func- 
tions 

Consider a theory with N colors and Nf massless flavors, with Nf small enough that the theory 
is asymptotically free. Label the flavors qf, qf. In an asymptotically free theory, we might hope 
that the short distance behavior of correlation functions can be systematically analyzed. The 
(gauge-invariant) Green's function 

Nf 

G{xi) = {l[qf{xf)qf{xf)) (1) 
/=l 

vanishes in perturbation theory, and receives a contribution from a single instanton. For sim- 
plicity, take one set of points, say xi,...,xa = x, and the other, xa+i, ■ ■ ■ XNf = 0, where 
Nf = 2A or Nf = 2A + 1, and call 

A Nf 
Oi{x) = llqfix)qf{x) 02{x)= J] (2) 
/=1 f=A+l 

We study the two point function 

A(x) = (Oi(a;)O2(0)). (3) 

We are interested in the behavior of A(x) for small x. If the correlation function is singular, 
we should be able to understand the singularities in terms of the operator product expansion. 
As for any pair of operators, the OPE of 01(0;), 02{x) takes the form: 

Oi (x)02 (0) « Yl ^"2 (^)^n (0) (4) 

The coefficient functions, in an asymptotically free theory like QCD, can be calculated in a short 
distance expansion. One expects that these include both perturbative and non-perturbative 
contributions. In the past, there has been much discussion of possible non-perturbative con- 
tributions to 0PE's[6, 5, 7]. In general, this is a subtle problem. If there are perturbative 



4 



contributions to a coefficient, any non-perturbative contribution whicli one miglit liope to iso- 
late will be exponentially small. In the case of A(a;), there are no perturbative contributions. 
One can ask: are any of the coefficient functions dominated by instantons? 

In perturbation theory, the theory respects an axial U{1) symmetry, as a result of which 
the lowest dimension operator allowed on the right hand side of eqn. 4 is simply OsNfi^^) = 
Oi (0)02(0), with coefficient C^^f{x). C^^f{x) starts out as a constant, and receives corrections 
in powers of as(x) (we will discuss these in section 5). More precisely, in perturbation theory, 
various operators of the same dimension as O^Nf appear on the right hand side, multiplied 
by powers of as{x)\ we will label these O3jy^(0), and their coefficients C^^^x). These will be 
considered more explicitly in section 5. Non-perturbatively, and in particular at the level of 
a single instanton, the unit operator can appear on the right hand side of this expression. In 
other words, 

Oi(x)O2(0) ~ C^^\x) + Cl^\x)0l^^{^) + higher dimension. (5) 

Assuming C^^"^ (x) is dominated by a well-defined instanton computation, we can determine its 
behavior on dimensional grounds (up to logarithmic variation associated with the perturbative 
anomalous dimensions of 0\ and O2): 

C(^) =A'^o^-3Nf+bo^ (6) 

where 60 - SNf = ^{N - Nf). So C(°)(a;) exhibits a power-law singularity at short distances 
if Nf > N, is non-singular if Nf < N, and possesses a logarithmic singularity for Nf = N. 
We will see that in the first case, the instanton computation is infrared finite, and there is a 
systematic expansion of the Green's function in powers of e 3^ and g^{x); in the second, 
the instanton computation is infrared divergent, and can be considered to correct the matrix 
elements of CsaTj ■ The Nf = N case requires more careful analysis. Correspondingly, there are 
three behaviors for the correlation function. For Nf > N, 

(Oi(x)O2(0)) « C{A''°x-^^f+''° + non singular) (7) 

where the coefficient C is calculable, but the non-singular terms are not. For Nf < N 

(Ci(x)C2(0)) « D{A^^f log(xA) + non singular) (8) 

where the coefficient, D, is not calculable. For Nf = N, the effects of instantons and other 
non-perturbative effects are comparable, and a more detailed analysis is required. This case, of 
course, is particularly interesting in real QCD, and will be the subject of section 5. 
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Note, most importantly, that in what we refer to as the non-calculable cases, the contribu- 
tions of the operator 01(0)02(0) fall off more slowly than the instanton contributions, so, for 
example, the contributions to this operator associated with dynamical breaking of the SU{Nf) 
chiral symmetry are more important than those calculated in the semiclassical approximation. 

Let's consider the instanton computation in more detail. In this section, we will examine 
the structure of the relevant integrals in order to determine their behavior with x; particular 
cases important for applications will be evaluated in more detail in section 7 and still more 
thoroughly in [8]. The fermion zero modes have the form: 

,i=p^ ^1 (9) 

where a is a two component spinor index and i is the gauge index. So the instanton contribution 
to A is 

A{x) = C / d^xodp- P (10) 

J [{x - xqY + p^Y^[xf^ + p^Y^'^ f ^1 

where C is a constant obtained from the non-zero modes in the instanton background, and xq 
and p are the translational and rotational collective coordinates, respectively. The integral over 
Xq may be performed by introducing Feynman parameters, yielding 

A(x) = C' j da[a'^-\l - af^^'^-^^-'^P^^^^^^^^ (11) 
J [x'^ai.l — a) -I- p^r-"*/ 



For large p, this behaves as 



I ^pf(^-^f). (12) 
P 



The integral converges for large p if Nf > N, exhibits a power law divergence for Nf < N, 
and diverges logarithmically for Nf = N. It is also free of ultraviolet divergences (which would 
appear at the endpoints of the a integration), and thus it behaves as guessed above: 

In the infrared divergent cases, the divergent part is identical to the (similarly ill-defined) 
instanton contribution to (Oi (0)02(0)). For Nf < N, the (cutoff) integral is non-singular for 
small X, corresponding to non-singular corrections to the coefficients of operators appearing in 
the OPE. For the case Nf = N, the expression also has a logarithmic singularity for small x, 
indicating the appearance of the unit operator in the OPE, with a coefficient function behaving 
as log(a;). It is necessary to define the operators appearing in these expressions at a scale M, 



^Note it is important that we divided the fields into two groups as we did; otherwise, there would be ultraviolet 
divergences associated with the definition of the local operators; this point is discussed further in section 8. 
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and this introduces a mass scale both into the matrix element and into the coefficient of the 
unit operator. We will discuss the scaling with M in more detail in section 5. Indeed, the 
question of whether we can compute the Green's function reliably is now a delicate one. The 
operators O^j^^^ possess anomalous dimensions. As a result, the coefficient functions behave as 
powers of logarithms. We see that the coefficient of the unit operator is proportional to a single 
power of a logarithm. So at best we can hope that C^^^ dominates by a power of a logarithm. 
We will explore this possibility in section 5. 



3 The Dilute Instanton Gas 



For the correlation functions we have been studying, there are two types of corrections we need 
to consider. First, there are perturbative corrections around the instanton. Examining the 
Feynman diagram expansion, these corrections are seen to be infrared finite, and so lead to 
corrections given by powers of as{x) to the coefficient of the unit operator. So we have, for the 
leading singular coefficient (displaying the result in a fashion which stresses the nature of the 
expansion). 

Here the first exponential reflects the anomalous dimensions of Oi and C>2- 

Now we ask whether there are corrections in powers of e "s^. Configurations containing 

— (2n+l)^4^ 

n + 1 instantons and n anti-instantons would be governed by e . As is well-known, in 

QCD there is not a controlled dilute instanton gas approximation, and we need to establish that 
such configurations do not some how swamp the calculable single-instanton contributions we 
have described above. The essential point of this section is that these dilute gas contributions 
to the Green's functions either represent corrections to the operator matrix elements, or non- 
singular corrections to the OPE coefficients. In other words, the complete expansion of A(x) 
is of the form of eqn. 13. 

A configuration with two instantons and one anti-instanton, widely separated, illustrates 
the main points. Specialize to the case N = 2, Nf = 3, and study the correlation function 
{Oi{x)O2{0)) ■ Each of the widely separated instantons or anti-instantons has six fcrmion zero 
modes. One can absorb these in various ways, leading to different types of contributions to the 
correlator: 
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1. "Disconnected Diagrams": Here all of the zero modes connected to one of the instantons 
are absorbed by the operators in the correlation function; the zero modes emerging from 
the other instanton and the anti-instanton are tied together (i.e. connected by propa- 
gators). The result is then the same as the one-instanton result, times an expression 
proportional to the space-time volume. This is a correction to the matrix clement of the 
unit operator (building the exponential e^^^^ in the Euclidean path integral). 

2. "Connected diagrams": Here there are contractions of fermion fields connecting the both 
instantons and the anti-instanton together. All such diagrams are infrared divergent, and 
do not exhibit singular behavior for small x. The leading divergences are x-independent, 
and correct the matrix elements of operators such as Oi (0)02(0). Subleading, non- 
singular x-dependent terms can correct the coefficient of the unit operator. 

Higher orders in the semiclassical expansion behave similarly. 

This argument establishes that widely separated instanton-anti-instanton configurations 
do not correct the singular part of the coefficient of the unit operator in the operator product 
expansion, and thus do not correct the singular behavior of A(x). Since the dilute gas expansion 
is not under systematic control in QCD, this can be viewed only as strong evidence that eqn. 
13 represents the complete expression for A(x), not a proof. It should be possible to test some 
of these questions in CP^ models [9]. 

4 The Large N Limit 

If N is large, and Nf fixed, the behavior of the correlation functions at short distances is readily 
understood using large N methods. Consider the case of two flavors, u and d. We are interested 
in the short distance behavior of: 



The matrix element factorizes in the large N limit, and so can be expressed in terms of {qq) 
(defined at the scale no). 



A{x) = (u(x)u(x)d(0)d(0)) 



(14) 



Prom the perturbative OPE, the right hand side includes a term: 




(15) 
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There will, in general, be additional operators on the right hand side, including the unit 
operator. We will comment on this possibility shortly. In any case, the expression above, in 
momentum space, behaves, at large momentum, as (it is more convenient to study the Fourier 
transform of V^A(a;), which is better behaved) 

V^A(x)|ft oc X powers of log(p). (16) 

Instantons give a finite contribution to the momentum space Green's functions, so these behave 
as 

5Ainst{x)^p^(-\ ; (17) 

for fixed k/p they are exponentially suppressed relative to (16) . If there are large N contri- 
butions to the unit operator in the OPE, (as one might reasonably suspect that there are), 
these may be even larger at large p. In any case, we have here a realization of the slogan that 
instanton contributions at large N are exponentially suppressed [4]. 



5 Nf = N 

We have seen that in the case A^^ = A^, the operators with dimension ?>Nf and the unit operator 
appearing in eqn. 5 in the OPE behave similarly with x, up to powers of logarithms. The best 
we might hope, then, is that the instanton wins (or loses) by a power of a logarithm. To assess 
this, we need to compute the anomalous dimensions of the various operators appearing in the 
OPE. This is a straightforward one loop computation. Consider, first, the case of SU{2) with 
two fiavors. Taking as the basis of dimension six operators 

Oi = uu dd O2 = uai^^u da^r,d (18) 

the matrix of anomalous dimensions is: 



where 



r has eigenvalues 



3/2 



A = -^-. (20) 



71 = ^2 = ^A (21) 
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Correspondingly, given that the /? function in this theory is 6, at small x, relative to the unit 
operator contribution, ln(a;), the contributions of the first operator behaves as (logx)^^/"^^, while 
those of the second behave as (logx)^/*^. So it is necessary to choose Green's functions carefully 
so as to obtain just the contribution of the second eigenoperator. 

In 5'C/(3), the matrix F is 3 x 3, and there are not simple analytic expressions for the 
eigenvalues, but one finds that, again, there is one operator which diverges more slowly at small 
X than the contribution of the unit operator. In order to isolate an instanton contribution, 
it is necessary to take a linear combination of operators in the Green's function, A(x), which 
projects on this. In practice, say in a lattice computation, this will be quite challenging, to say 
the least. 



6 Application: Correlation Functions for e^e annihilation 

Given our experience with chirality violating Green's functions, it is interesting to ask whether 
we can isolate instanton effects in physical processes. This is a question which has been inves- 
tigated in the past (sec, for example, [5, 6, 7]). Consider, for example, e+e~ annihilation (well 
below the threshold). The hadronic total cross section is proportional to 

T{J>^{x)r{y)), (22) 

Fourier transformed to momentum q. The challenge is to separate out instanton effects from 
other perturbafAve and non-perturbative effects. The framework wc have developed here is useful 
for this. Before performing any calculations, wc note the following. In perturbation theory, the 
unit operator appears in the operator product expansion, with a coefficient function which 
behaves as 1/x^. Fourier transformed, this gives the familiar form: 

U^"'{q) = {q^gi''' - q>'q'')n{q^) (23) 

with the perturbative contributions to 11 behaving as powers of logarithms. There is no reason 
to think that there aren't power law corrections at large q^, proportional to powers of A^/g^. 

Let's consider instanton contributions. More precisely, we need to consider dilute gas 
contributions with net instanton number zero. In general, these will be ir divergent in coordinate 
space, unless 6o < 3/2, and only for 6o = are they as singular at short distances as the 
perturbative contributions. The leading IR divergences can be thought of as corrections to the 
coefficient function of the unit operator (and higher dimension operators). In momentum space. 
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they will fall rapidly with momentum, and thus cannot be distinguished from other types of 
non-perturbative effects (along the lines of our large N argument above). 

Taking the lessons from our earlier discussion, we can establish the result by dimensional 
analysis. The p integral, in momentum space, is always convergent [6], so the result of an 
instanton-anti-instanton contribution is necessarily: 

n(g2) = A2^'Og-2bo_ (24) 

To be larger than the perturbative contributions, we require, again, 60 < Oj i-e. we need 
be outside the regime where the theory is asymptotically free. If one examines the actual 
amplitudes, one finds expressions in agreement with that above (the coordinate space expression 
is infrared finite for <2). Of course, in large we have an exponential suppression. 

So, at least in this context, it does not seem possible to isolate instanton contributions to a 
physical, high momentum cross section. Similar remarks apply to deep inelastic scattering and 
processes in heavy quark physics. To find computable physical quantities requires, at the very 
least, a theory in which chirality operating operators arc important observablcs. We do not sec 
an obvious application to realistic scattering experiments. However, as we describe in the next 
section, there are potentially important applications of these ideas in lattice gauge theory. 

7 Applications: Lattice Gauge Theory in the Chiral Limit 

A potential application of these ideas lies in tests of lattice gauge theory simulations. What we 
have just learned is that, for suitable numbers of flavors and colors, maximally chirality violating 
correlation functions can be computed systematically in a semiclassical approximation. Given 
that the study of the chiral limit of lattice gauge theories is subtle, while, at the same time, 
current analyses probe quite small quark masses [10], comparison of "experiment" and "theory" 
could be useful. 

The realistic case, for which a significant amount of data exists (i.e. for which there already 
exists a large ensemble of gauge configurations), would he Nf = N = 3. We have seen that 
this is the borderline situation for the semiclassical analysis; at the same time, the large N 
analysis is not reliable. At most, the instanton computation dominates by a small power of 
a logarithm. Still, it would be of interest to examine the data to see if the Green's function 
exhibits roughly logarithmic variation with a;, of a suitable order of magnitude. The simplest 
case to study where a reliable computation is possible is the case N = 2,Nf = 3. While not 
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nature, this should be accessible to detailed numerical study. Here 

A16/3 

(Oi(x)O2(0)) = C-^{1 + 0{as{x))). (25) 

where the numerical coefficient, C, is presented below in a particular scheme. The question, 
then, is whether, with current lattice technology, such a computation is feasible; as we now 
explain, we believe that, while such a calculation is challenging, given the lattice spacings and 
quark masses currently being explored, the answer may well be yes. A detailed analysis of the 
issues will appear in reference[8]; here we content ourselves with some rough estimates. 

In practice, one would compute the correlation function, dividing out by the vacuum to 
vacuum amplitude ((I)). Current lattice computations achieve quite small lattice spacings and 
quark masses.. For non-zero quark mass, perturbing in the quark masses, there are more 
singular contributions to the correlation function than those due to instantons. (The relevant 
parameters are m^/A and l/(Ax), where A is the QCD scale). 

First, we estimate more carefully the size of instanton effects. The needed functional 
determinants can be obtained from 't Hooft's original paper [1]. Here we will content ourselves 
with a brief summary of special cases in the MS scheme; more general results, including lattice 
regulators, will appear in ref. [8]. For the case Nf = 3,N = 2, we obtain: 

A(x) « 9 X lO^n^'^/^g-^e^x'^^/^ (26) 
while the logarithmically singular term, for Nf = 3,N = 3 is: 

A(a;) = -1.0 X 10V^5"^^e^ log(M x). (27) 



To estimate the effects on the lattice size, we will take the smallest lattice spacing and 
quark masses considered by the MILC collaboration[10]; roughly a = (4 GeV)~^, and m^, Ri 
10 — 20 MeV. We will consider measuring the correlation function at a; = (1.5 GeV)~^ (by these 
scales, we mean scales such that the SU (2) gauge coupling has the value corresponding to the 
observed SU{3) coupling at that scale). The leading terms are independent of x and quite 
sensitive to the cutoff, a; roughly they give a contribution, in the case of three flavors: 

A=^^^. (28 

Using the expression for A, eqn. 26, this is overwhelmingly larger than the instanton contri- 
bution. This problem can be avoided by choosing operators which are odd under parity (the 
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parity symmetry in the case of staggered fermions is discussed, for example, in [11]). In that 
case, the leading contribution arises from a four loop diagram. This behaves as (a detailed 
computation will be presented in [8]): 



and should be much smaller than the instanton contribution. 

One then can ask whether the instanton contribution is large enough that it is realistic 
to hope to extract it. We believe the answer is yes, but that the computation is likely to be 
challenging. A yardstick (fermistick?) for comparison is 



with some contraction of the indices on Sp. This is roughly comparable in size to A for 
= 1.5 GeV. So one might hope that the numerical computation will not be so noisy as to 
mask the instanton contribution. Similar statements apply to the SU (3) computation (which, 
again, is subject to greater uncertainties, but is clearly of interest). Examining somewhat 
similar numerical computations reinforces our optimism that this is a challenging, but tractable, 
computation[12]. 

8 Other OPE's, and the Definition of the Composite Operators 

In section 2, we studied Green's functions for operators where we split the 2Nf fermion fields 
more or less evenly. If we do not do this, in general, for Nj > + 1, we encounter ultraviolet 
divergences. These can, again, be usefully understood in the language of the operator product 
expansion; if we define the local operators by point splitting, then, under some circumstances, 
the instanton contribution to the point split operator diverges as the splitting tends to zero. 

Consider the N = 2,Nf = 3 case. Here we would expect that the OPE of the operators 
uu{x)dd{0) includes s*s* (in two component notation). The question is whether the OPE 
coefficient is singular. It is easy to check that the three point function. 



in the limit |xp ^ is non-singular. So the corresponding OPE coefficient, while non-zero, 
is also non-singular as x ^ 0. More precisely, the OPE of pairs qiqi includes terms: 




(29) 




(30) 



{u{x)u{x)d{0)d{0)s{y)s{y)) 



(31) 



u{x)u{x)d{0)d{0) 



u(0)u(0)d(0)d(0) + f{x)s*{0)s*{0) + ... 



(32) 
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The statement that the Green's function of eqn. 31 is non-singular for smaU x is the statement 
that the coefficient function, f{x), is non-singular for small x. The leading term in the Green's 
function, indeed, is precisely the matrix element of the first operator in eqn. 32. 

As in the case of the unit operator we encountered earlier, if we increase the number of 
flavors, for fixed N, the analog of f{x) becomes singular. For example, for N = 2,Nf = 4, one 
can study the OPE: 

u{x)u{x)d{x)d{x)s{0)s{0) ~ n(0)u(0)(i(0)(i(0)s(0).s(0) + f {x)t* {0)t* (0) + ... (33) 

where now f{x)A~^^/^ ^ x~^^^. Again, the coefficient function can be calculated using instan- 
tons. 

So, for example, for four flavors in SU{2), if we study the Green's function: 

{u{x)u{x)d{x)d{x)s{x + e)s{x + e)f(0)t(0)) (34) 

and simply set e = 0, we obtain an ultraviolet divergence, associated with the singularity of the 
operator product expansion discussed above. 

9 Conclusions and Directions for Further Research 

We have seen that scmiclassical methods are sometimes applicable to QCD, including cases with 
realistic numbers of flavors and colors; we have also seen that for many questions, as anticipated 
long ago, instanton methods are not useful[4]. There are Green's functions which can be 
computed systematically using instanton methods if Nf > N. The case Nf = N is borderline, 
in which, at most, the instanton dominates at short distances by a power of a logarithm. For 
iV S> iV, we have seen an explicit realization of the expectation that instanton effects are 
exponentially at small N. On the other hand, the most interesting case of Nf = N = 3 exhibits 
no such suppression. 

An interesting arena in which to analyze all of these questions is provided by the CP^ 
models. For such theories, an exact large N solution is available, and it should be possible 
to test a number of the ideas developed here. It would, however, be especially interesting to 
attempt to study analogous Green's functions at finite N. These issues will be discussed in a 
subsequent publication [9]. 

The most immediate application of these results is likely to be in lattice gauge theory. 
Given the state of the art of current lattice computations, which probe quite short distances 
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and small quark masses, one might hope to extract some of the Green's functions discussed 
here, comparing not only their qualitative behavior but their precise numerical values. Results 
for the relevant Green's functions, including the effects of finite quark mass, will be presented 
in [8]. Successful agreement would provide strong evidence for the accuracy of such lattice 
computations, and make important assertions, such as the fact that niu 7^ 0, all the more 
convincing. 
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